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SYMBOLS 


A  nozzle  area 

a  speed  of  sound 

Cp  particle  drag  coefficient 

d  diameter 

h  enthal py 

M  Mach  number 

m  mass 

p  pressure 

S  reference  area 

s  distance 

t  time 

V  velocity 

v  reduced  velocity,  V/Vm 

x  non-dimensional  distance,  s/A 

Y  ratio  of  specific  heats 

A  characteristic  length  defined  by  Equation  (7) 
p  density 

Subscripts  and  Superscripts 

f  nozzle  exit 

m  gas  stream  maximum 

p  particle 

s  gas  stream 

t  gas  stagnation  conditions 

*  sonic  flow  conditions 
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Section  1 


INTRODUCTION 


The  rational  design  of  a  hypersonic  vehicle  requires  the  ability  to  predict 
the  behavior  of  the  vehicle  operating  within  the  environment  to  which  it  is 
exposed.  This  environment  may  include  high  speed  flight  through  an  atmosphere 
containing  rain,  ice  and/or  dust  particles  which  results  in  both  particle  erosion 
as  well  as  ablation  of  the  heat-protection  system.  Neglecting  the  erosion  due 
to  particles  can  result  in  destruction  of  the  vehicle  with  resultant  failure 
of  the  mission.  Computer  programs  of  various  degrees  of  sophistication  have 
been  developed  at  MDAC,  and  by  other  investigators,  to  predict  the  behavior  of 
thermostructural  systems  in  such  environments.  These  programs  require  empirical 
erosion  coefficients  which  have  been  determined  from  simulated  ablation-erosion 
experiments. 

Existing  hypervelocity  erosion  facilities  do  not  simulate  the  real  flight 
environment  completely  and  particularly  for  the  advanced  vehicles.  The  short¬ 
comings  may  consist  of  one  or  more  of  the  following  characteristics:  too  low 
velocities,  too  short  exposure  times,  lack  of  simultaneous  heating  and  erosion 
simulation,  inadequate  means  of  dispersing  the  particles,  inadequate  size  of 
model  or  particles  or  both.  Some  of  the  methods  used  to  obtain  erosion  data 
include  the  following:  (1)  accelerating  a  small  group  of  particles  by  means 
of  a  shock  tunnel  ("hot-shot")  or  a  light-gas  gun  and  impinging  these  particles 
upon  a  specimen  which  has  been  preheated  by  a  plasma  arc  or  by  a  radiant  heater; 
(2)  firing  the  specimen  down  a  ballistics  range  (maximum  test  length  less  than 
1,000  ft)  and  recording  data  in  flight;  (3)  rocket  sled  testing  which  is  expen¬ 
sive,  Is  limited  to  velocities  below  flight  velocities,  and  has  not  been  developed 
to  the  point  of  obtaining  erosion  data  with  small  dust  particles;  (4)  hypersonic 


wind  tunnels  Into  which  dust  particles  are  Injected  provide  more  realistic 
testing  times  and  multiple  Impingement  but  low  velocities,  low  heating,  and 
low  shear.  All  these  methods  provide  widely  varied,  and  limited  levels  of  erosion 
simulation.  A  definite  need  exists  for  more  realistic  simulation  of  the  flight 
environment  which  can  result  in  combined  ablation  and  particle  erosion. 

Of  the  current  methods  for  conducting  erosion  tests  only  the  wind  tunnels  can 
provide  testing  times  of  sufficient  duration  to  match  flight  times  (of  the 
order  of  10  seconds).  The  characteristics  of  wind  tunnels  which  result  in  low 
velocities  and  low  heating  are  the  relatively  low  values  of  stagnation  enthalpy 
and  stagnation  pressure  associated  with  such  facilities.  These  limitations 
may  be  overcome  by  utilizing  another  facility  as  the  energy  source  for  a 
particle  accelerator,  namely,  the  plasma  arc  which  has  been  used  extensively 
for  ablation  testing.  In  particular,  if  we  Inject  the  dust  particles  at  or  near 
the  throat  of  a  high-enthalpy,  high-pressure  plasma  arc,  then  utilize  a  high 
Mach  number  nozzle  to  accelerate  the  gas  stream  to  high  velocities,  the  gas 
stream  will  accelerate  the  dust  particles  to  velocities  that  can  approach  the 
stream  velocity.  A  facility  of  this  type  would  provide  a  reasonable  testing 
time  together  with  higher  particle  velocities  and  larger  heating  rates  than 
current  wind  tunnels. 

Before  embarking  on  the  design  and  hardware  effort  of  modifications  to  existing 
arc  facilities,  it  was  felt  that  an  analytical  Investigation  was  mandatory  to 
explore  the  phenomena  involved  in  the  energy  transfer,  from  the  internal  energy 
of  the  gas,  to  the  kinetic  energy  of  the  particles.  A  primary  objective  was 
to  Investigate  methods  of  tailoring  the  gas  flow  In  order  to  obtain  efficient 
acceleration  of  the  particles,  in  particular,  to  determine  methods  for  optimizing 


such  gasdynamic  particle  accelerators.  The  results  of  this  investigation  should 
provide  the  methods  for  an  assessment  of  the  capabilities  of  existing  wind 
tunnels  and  arc  jets  as  particle  accelerators  and  would  provide  the  means  for  a 
rational  approach  to  the  modification  of  existing  facilities  as  well  as  the 
design  of  new  facilities  for  conducting  steady-state  dust  erosion  experiments  and 
combined  ablation-erosion  experiments.  Moreover,  the  availability  of  an  optimum 
solution,  even  if  somewhat  idealized,  is  useful  for  many  reasons.  Some  of 
these  include:  a  goal  to  which  a  designer  can  strive;  a  basis  for  evaluating 
an  existing  design;  a  focus  upon  the  limitations  that  nature  imposes  upon  a 
problem;  and,  more  often  than  not,  an  indication  of  the  modifications  that 
would  improve  the  performance  of  an  existing  design. 
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Section  2 
ANALYSIS 


2.1  Equations  of  Motion 

The  basic  equation  of  motion  for  a  particle  in  a  gas  stream,  such  as  a  wind 
tunnel,  arc- jet,  or  nozzle,  may  be  written  in  the  form 


C0  0,  Sp  <vs  -  V‘ 


"p 

m  — t 
p  dt 


(1) 


Equation  (1)  may  also  be  written,  by  noting  that  dVp/dt  *  (dVp/ds)(ds/dt)  « 
Vp  dVp/ds  where  s  Is  the  arc  length  along  the  trajectory  of  the  particle. 


cp  Ps  SP  (V,  -  V 


dV 

I)  Vn  — E 

P  P  ds 


(2) 


The  general  one-dlmenslonal  problem  of  particle  motion  in  a  gas  stream  corresponds 
to  specifying  the  gas  density,  p$(s),  and  the  stream  velocity,  V$(s),  as  arbi¬ 
trary  functions  of  distance  along  the  gas  stream  and  then  Integrating  equation 
(2)  to  obtain  the  particle  velocity,  Vp,  as  a  function  of  the  distance,  s,  or 
of  the  time,  t,  (In  equation  [1]). 


For  a  perfect  gas,  with  Isentropic  flow,  the  gas  density  may  be  related  to  the 
stream  velocity,  V$,  and  the  stagnation  density,  pt,  through  the  relation 
(reference  1) 


(3) 
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Substituting  equation  (3)  into  equation  (2)  and  dividing  by  gives 


) 


(4) 


Note  that  for  any  gas  stream  (with  specified  stagnation  conditions)  there  exists 
a  maximum  velocity,  Vm,  that  is  assocaited  with  the  given  stream.  The  velocity 
corresponds  to  the  velocity  attainable  by  exapnding  the  gas  to  zero  pressure  and 
is  related  to  the  stagnation  speed  of  sound,  a^,  by  the  relation 


(5) 


Substituting  equation  (5)  into  equation  (4)  gives 

Lz-pji  uj  1  \vm  vm;  ds 


(6) 


Substituting  the  nondimensionalized  velocities  vs  and  vp,  and  noting  that  the 
quanti ty 


X  - 


2  m_ 


CDSppt 


(7) 


has  the  dimension  of  length,  we  can  nondimensionalize  the  distance  s  by  dividing 
by  X,  so  that  x  *  s/X  and  equation  (6)  becomes  (when  X  is  a  constant) 


(8) 


The  problem  of  evaluating  the  effectiveness  of  a  gasdynamic  stream  for  accel¬ 
erating  particles  is  Imbedded  within  the  solutions  of  equation  (8).  The 
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range  of  variables  of  interest  consists  of  nondimensional  velocities  bounded 
by  0  <  Vp  £  vs  <  1.  The  direct  problem  consists  of  obtaining  solutions  to 
equation  (8)  for  specified  variations  of  the  nondimensional  gas  velocity  v$(x) 
as  a  function  of  the  distance.  It  should  be  noted  that  equation  (8)  can  be 
modified  to  include  the  case  of  gasdynamic  decelerators  (vp  >  v$)  by  multiplying 
the  right  hand  side  by  minus  one.  It  should  also  be  noted  that  rather  large 
ranges  of  physical  situations  are  included  in  the  apparently  simple  format  of 
equation  (8).  These  physical  conditions  include  the  following:  arbitrary 
particle  size  (e.g.,  from  less  than  1  micron  in  diameter  to  well  over  1  centimeter); 
the  complete  range  of  stagnation  conditions,  such  as,  pressure,  density,  and 
pressure  (within  the  range  of  a  perfect  gas);  and  the  full  range  of  stream  Mach 
number  from  zero  to  infinity. 

A  preliminary  assessment  of  the  character  of  equation  (8)  which  governs  the 
motion  of  particles  in  an  accelerating  gas  stream  indicates  that  its  apparent 
simplicity  is  deceiving.  The  equation,  although  of  first  order,  is  nonlinear 
and  with  variable  coefficients.  The  analytical  solution  of  such  equations  is 
not,  in  general,  an  easy  task.  The  particular  equation  with  which  we  are  con¬ 
fronted  here  can  be  shown  to  be  of  the  Riccati  type.  Such  equations,  in  the 
general  case,  are  not  reducible  to  quadratures.  One  approach  to  the  analytical 
solution  of  the  basic  equation  of  particle  motion  in  a  gas  stream  would  be  to 
make  judicious  approximations  in  order  to  render  the  equation  tractable  while 
retaining  the  essence  of  the  problem.  This  method  was  used  in  reference  2  and 
reasonable  results  appear  to  have  been  obtained.  Another  approach  is  to  seek 
exact  solutions  to  special  problems.  In  this  approach  two  categories  of  special 
problems  are  available.  The  first  is  to  seek  solutions  to  indirect  problems  in 
which  the  particle  motions  (vD(x))  are  specified  and  then  to  solve  for  the  stream 
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velocity  ( vg ( x ) )  required  to  obtain  the  given  motion.  This  approach  circumvents 
the  solution  of  the  nonlinear  differential  equation;  however,  the  success  of 
this  method  depends  upon  the  astuteness  of  the  specified  motion  and  its  relation 
to  physical  reality.  The  more  interesting  approach  is  to  seek  optimization 
solutions.  The  present  investigation,  therefore,  focussed  upon  the  flowing 
optimization  problem:  For  every  particle  velocity,  vp,  find  the  value  of  stream 
velocity,  v$,  that  maximizes  the  acceleration  of  the  particle.  Such  a  solution 
gives  the  shortest  distance  to  accelerate  the  particles  from  one  value  of 
velocity  to  another  value. 

2.2  Optimization 

In  order  to  determine  the  stream  velocity,  v$,  that  maximizes  the  acceleration 
for  every  given  particle  velocity,  the  left  hand  side  of  equation  (8)  is 
differentiated  with  respect  to  vs>  and  the  derivative  is  set  equal  to  zero.  The 
resulting  equation  is 

n-vs2)  T-1  (vs-vp)Z  +  2n-vs2jY"1  (Vvp>  «  0  (9) 

One  solution  to  equation  (9)  is  v$  ■  vp,  the  uninteresting  case  for  which  the 
particle  acceleration  vanishes.  The  other  solution  is  readily  obtained  from 
equation  (9)  by  solving  for  vp,  thus, 

vp  =  YVs  "  0°) 


Equation  (10)  provides  the  desired  relation  between  the  dimensionless  stream 
and  particle  velocities  for  which  the  particle  acceleration  is  a  maximum.  The 


dv 

associated  maximum  nondi  mens  tonal  particle  accleratlon,  vp“g7  mx*  Is 
obtained  by  substituting  equation  (10)  Into  equation  (8)  to  give 


fv 

\  p  dx/max 


(11) 


Although  equation  (10)  represents  a  significant  step  In  the  optimization  problem 
by  relating  the  two  velocities,  there  reamlns  to  relate  these  velocities  to  some 
geometric  parameter  such  as  the  nondi mens ional  distance,  x.  In  the  direction  of 
the  flow;  i.e.,  the  functions  v$(x)  and  vp(x)  are  to  be  determined.  One  approach 
Is  to  solve  equation  (10)  for  v$,  substitute  the  result  Into  equation  (8)  to 
obtain  a  differential  equation  for  vp,  and  then  solve  the  differential  equation. 
This  procedure  results  In  Integrals  that  are  difficult  to  evaluate.  A  more 
fruitful  approach  was  found  to  be  to  transform  equation  (10)  Into  a  differential 
equation  in  v$  in  the  following  manner.  Differentiate  equation  (10)  with  respect 
to  x  to  obtain 


dv„  dv„ 


oa) 


Then,  substitute  equations  (10)  and  (12)  Into  equation  (8)  to  obtain,  after 
some  algebraic  manipulation,  the  following  differential  equation 


(13) 


This  differential  equation  can  be  solved  for  x  by  separation  of  variables  to 
give 
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Thus,  the  problem  of  determining  the  variation  of  dimensionless  stream  velocity 
with  dimensionless  distance  along  the  flow  direction  to  obtain  maximum  particle 
acceleration  has  been  reduced  to  a  quadrature.  The  associated  particle 
velocities  may  then  be  obtained  from  equation  (10). 

The  Integrals  in  equation  (14)  have  been  evaluated  using  the  tables  of  references 
(3)  and  (4)  for  three  constant  values  of  the  ratio  of  specific  heats,  y;  namely, 
5/4,  7/5,  and  5/3.  These  values  of  y  are  representative,  respectively,  of  air 
in  the  nozzle  of  a  high  temperature  air  arc,  a  perfect  diatomic  gas,  and  a 
perfect  monatomic  gas.  After  considerable  algebraic  manipulation,  the  expressions 
for  x  as  a  function  of  v$  were  reducible  to  the  following  forms: 

for  y  »  5/4: 
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Section  3 


RESULTS  AND  DISCUSSION 


The  relation  between  reduced  particle  velocity.  vp,  and  reduced  stream  velocity, 
v$,  for  maximum  particle  acceleration  is  shown  in  Figure  1  for  several  values 
of  the  ratio  of  specific  heats,  y.  Although  physically  unrealistic,  the  case 
of  y  *  1  is  included  to  delineate  a  limiting  case  which  corresponds  to  zero 
particle  acceleration  (since  v$  ■  vp).  The  value  of  v$  when  vp  *  0  is  of  some 
Interest  since  it  corresponds  to  the  stream  velocity  that  gives  maximum  accel¬ 
eration  to  the  particle  when  the  particle  is  Injected  Into  the  stream  with  zero 
initial  velocity.  In  this  particular  case  the  dynamic  pressure  acting  upon 
the  particle  is  Identically  equal  to  the  dynamic  pressure  of  the  gas  stream. 

This  value  of  v$  Is  readily  obtained  from  equation  (10)  by  setting  vp  *  0  to 
give 


v 


s 


when  v  ■  0 
P 


(18) 


Since  the  Mach  number,  Ms,  of  the  gas  stream  is  related  to  the  reduced  stream 
velocity,  vs,  by 


(19) 


It  is  easily  shown  that 


Ms  ■  JT  when  vp  ■  0  (20) 

for  maximum  particle  acceleration.  Moreover,  this  result  Is  valid  for  all 
values  of  the  ratio  of  specific  heats,  y.  Thus,  in  order  to  maximize  particle 


acceleration,  the  appropriate  location  for  injecting  particles  with  zero  initial 
velocity  Into  a  nozzle  is  the  point  at  which  the  gas  stream  Mach  number  is 
equal  to  *7. 


dv 


LVP 


“at. 


max 


,  with 


The  variation  of  reduced  maximum  particle  acceleration, 
reduced  stream  velocity  is  presented  In  Figure  2  for  three  values  of  y.  Also 
Included  In  Figure  2  is  a  dashed  curve  for  the  variation  of  maximum  particle 
acceleration  when  the  particle  velocity  Is  zero.  The  expression  for  this  curve 
can  be  derived  by  substituting  equation  (18)  Into  equation  (11)  to  obtain 


[VP  -a&]«x - 


when  v_ 


(21) 


Equation  (21)  gives  the  largest  values  of  reduced  particle  acceleration  for  the 
problem  considered  In  the  present  Investigation;  namely,  the  acceleration  of 
particles  with  zero  or  finite  velocity  in  the  direction  of  a  flowing  gas.  It 
Is  possible  to  have  larger  magnitudes  of  particle  acceleration  by  injecting 
the  particles  upstream  Into  the  flowing  gas  (vp<0);  the  deceleration  of  particles 
Is  not,  however,  the  concern  of  the  present  Investigation.  It  will  be  noted 
from  Figure  2  that  the  reduced  acceleration  Increases  with  Increasing  y.  Also 
from  Figure  2  It  Is  clear  that  the  maximum  particle  acceleration  decreases 
rapidly  as  the  reduced  stream  velocity  approaches  the  limiting  gas  velocity 
(v^l.O).  In  view  of  these  low  values  of  maximum  particle  acceleration  it 
would  be  expected  that  nozzles  designed  as  particle  accelerators  will  become 
much  longer  as  one  requires  the  particle  velocity  to  approach  more  closely  the 
limiting  gas  velocity. 


Another  Illustration  of  the  variation  of  reduced  particle  acceleration  with 
reduced  velocities  (particle  velocity  and  stream  velocity)  Is  shown  In 


Figure  3  where  lines  of  constant  (v  -g£)  are  plotted  on  the  (vp,  vs)  -  plane 
for  y  «  7/5.  Also  Included  in  Figure  3  is  the  optimum  variation  of  vp  and  v$ 
(repeated  from  Figure  1).  It  will  be  noted  that  the  acceleration  is  flat  in 
the  neighborhood  of  the  optimum  curve.  In  addition  the  acceleration  vanishes 
when  vs  *  1  and  when  vp  ■  v$.  Thus,  a  nozzle  design  that  deviates  significantly 
from  the  optimum  curve,  either  to  the  right  or  left  of  the  optimum,  will  be 
significantly  less  effective  in  accelerating  particles.  In  the  region  above 
and  to  the  left  of  the  optimum  curve  the  particle  slip  velocity  (vs"vp)  1*  too 
low  although  the  stream  density  retains  a  slgnificent  value.  In  the  region 
below  and  to  the  right  of  the  optimum  curve  the  nozzle  is  overexpanded  and  the 
density  is  too  low,  although  the  slip  velocity  may  remain  large.  In  either 
case,  the  dynamic  pressure  acting  on  the  particle  is  less  than  the  maximum. 

Inasmuch  as  the  curves  for  constant  reduced  acceleration  are  flat  in  the 
neighborhood  of  the  optimum,  minor  deviations  from  the  optimum  may  not  result 
In  significant  reductions  in  particle  acceleration.  Plots  of  the  form  of 
Figure  3,  l.e.,  the  (vp,  v$)  -  plane,  may  be  used  as  a  guide  for  modifying  an 
existing  nozzle  design  to  improve  its  ability  to  accelerate  particles.  In 
particular,  if  one  plots  the  performance  of  a  given  design  as  a  curve  on  the 
(vp»  vs)  ■  plane  and  compares  It  with  the  optimum  curve,  the  qualitative 
changes  In  nozzle  design  become  apparent:  expand  more  rapidly  those  sections 
of  the  nozzle  that  lie  to  the  left  of  the  optimum  curve  and  reduce  the  expan¬ 
sion  ratio  for  the  parts  of  the  nozzle  represented  by  a  curve  to  the  right  of 
the  optimum. 

The  variation  of  the  reduced  velocities,  v$  and  vp,  with  dimensionless  distance, 
x,  along  the  flow  direction  to  obtain  maximum  particle  acceleration  Is  presented  In 
Figure  4  for  three  values  of  the  ratio  of  specific  heats,  y.  equal  to  5/4,  7/5, 


and  5/3.  The  semi-log  format  of  Figure  4  is  used  to  spread  out  the  variations 
for  small  values  of  x  where  the  accelerations  are  large  with  resulting  rapid 
changes  in  velocity.  It  will  be  noted  from  Figure  4  that  large  Increases  in  nozzle 
length  are  required  for  modest  increases  In  particle  velocity  when  the  reduced 
velocities  approach  unity.  The  larger  values  of  reduced  velocities  with 
increasing  y  that  were  noted  previously  in  Figure  2  can  be  seen  to  result  in 
significantly  smaller  values  of  x  for  the  same  reduced  particle  velocity.  The 
optimum  streamwise  variation  of  the  dimensionless  parameters  v$  and  vp,  together 
with  gas  stream  Mach  number,  Ms,  and  nozzle  area  ratio,  A/A*,  are  also  shown  In 
Figure  5  using  linear  coordinates  for  the  case  of  y  •  7/5. 

Some  typical  variations  of  nozzle  parameters  (d  /d  *,  v  ,  v  ,  and  M  )  with 
the  dimensional  distance,  s  (In  Inches),  downstream  from  the  particle  Injection 
point  are  presented  In  Figure  6  through  9  for  high-temperature  air  nozzles 
that  maximize  the  particle  acceleration.  The  stagnation  conditions  of  the  air 
were  taken  to  be  :  68  atmospheres  pressure,  591 3#R  temperature,  and  0.458  lbm/ 
ft3  density.  The  value  of  the  ratio  of  specific  heats,  y,  was  taken  to  be  equal 
to  5/4,  which  Is  representative  of  high- temperature  air  In  an  arc.  The  particles 
were  magnesium  oxide  with  a  density,  pp,  equal  to  223.5  lbm/ft3  and  particle 
diameters,  dp,  equal  to  5,  50,  and  500  microns.  The  particle  drag  coefficient, 

Cq,  was  constant  and  equal  to  1.0  for  most  of  the  curves;  a  value  of  Cp  ■  2.0 
was  also  used  with  the  50  micron  case  to  Indicate  the  effect  of  changing  particle 
drag  coefficient.  The  results  of  Figures  6  through  9  show  that  there  exists  an 
optimum  nozzle  for  each  particle  size  and  that  as  particle  size  Increases  the 
length  of  the  optimum  nozzle  to  achieve  a  given  particle  velocity  also  Increases. 
This  result  Is  a  consequence  of  the  fact  that  the  ratio  of  Inertial  to  drag  force 
Is  proportional  to  \  which,  In  turn.  Is  proportional  to  particle  diameter.  Thus, 
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It  Is  more  difficult  to  accelerate  large  particles  by  a  gas  stream  than  it  is  to 
accelerate  the  smaller  ones.  In  general,  changing  those  parameters  which  decrease  X 
will  result  In  shorter  optimum  nozzles  for  a  given  particle  velocity  or  for  obtaining 
higher  particle  velocities  from  a  given  nozzle  length.  From  the  definition  of  X 
it  can  be  shown  that  decreasing  particle  size  and  particle  density  as  well  as  in¬ 
creasing  particle  drag  coefficient  and  stagnation  density  of  the  gas  will  decrease  X. 

A  meaningful  comparison  can  be  made  with  the  results  of  the  present  investigation 
and  some  recently  published  calculations  of  Reference  2  for  the  case  of  a  conical 
nozzle.  This  comparison  is  shown  in  Figure  10  where  the  particle  velocity, 

Vp  ,  at  the  exit  of  the  nozzle  is  plotted  as  a  function  of  particle  diameter  dp. 

The  stagnation  conditions  of  the  flow  were:  75  atmospheres  pressure,  2000  Btu/lbm 
enthalphy.  The  maximum  stream  velocity,  Vm,  for  this  value  of  stagnation  enthalpy 
was  found  to  be  10,000  feet  per  second  when  using  the  energy  equation  for  air 
given  in  Reference  5  as 

V$  *  223.6  (ht  -  h)1/2  (22) 

where  the  velocity  Is  given  in  feet  per  second  and  the  enthalpies  are  in  Btu/lbm. 
Setting  h  equal  to  zero  In  equation  (22)  gives 

«  223.6  (ht)  1/2  (23) 

The  nozzle  length  was  17.25  feet  and  both  sets  of  calculations  assumed  a  constant 
drag  coefficient  CD  *  1.0.  The  differences  in  ratio  of  specific  heats  for  the 
two  sets  of  calculations  [y  *  1.22  in  Reference  [2]  and  y  •  1.25  for  the 
optimum  nozzle)  are  deemed  insignificant.  It  will  be  noted  from  Figure  10  that 
whereas  for  small  particle  sizes  only  minor  Improvements  In  particle  velocity 
are  possible  (approximately  800  feet  per  second  increase  over  6000  feet  per 
second  for  50  micron  diameter  particles),  significant  increases  are  possible 
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for  the  larger  particle  sizes  (4040  feet  per  second  compared  to  approximately 
2400  feet  per  second  for  1500  micron  diameter  particles).  It  should  be  noted 
that  the  stagnation  conditions  chosen  for  this  calculation  do  not  represent 
the  current  state-of-the-art  In  plasma  arc  facilities. 

Another  method  for  Increasing  the  particle  velocities  Is  to  Increase  the  value 
of  the  maximum  stream  velocity,  V^,  of  the  gas.  From  equation  (23)  It  Is  clear 
that  the  stagnation  enthalpy  of  the  gas  must  be  Increased.  This  effect  is 
illustrated  by  the  upper  curve  of  Figure  10  which  shows  the  substantial  increases 
of  particle  velocities  in  an  air  arc  facility  with  a  stagnation  enthalpy  of 
6000  Btu/lbm  (at  190  atmospheres  stagnation  pressure  in  order  to  keep  X 
essentially  Invariant  between  the  two  sets  of  calculations).  For  this  case 
the  calculations  give  (for  CD  -  1.0)  particle  velocities  of  11,800  and  7,000 
feet  per  second  for  particle  diameters  of  50  and  1,500  microns,  respectively. 

Enthalpies  at  this  pressure  level  appear  to  be  an  attainable  goal  for  the  near 
future  (Reference  6).  Thus,  a  plasma  arc  with  an  appropriately  optimized 
nozzle  appears  to  offer  the  possibility  of  significantly  improved  simulation 
of  the  high  heating  rates  simultaneously  with  the  high  particle  velocities  and 
longer  testing  times  required  for  advanced  vehicles. 

\ 
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Section  4 


CONCLUSIONS 

The  problem  of  determining  the  streamwise  variation  of  gas  stream  velocity  in 
a  nozzle  such  that  the  acceleration  of  particles  within  the  gas  stream  is 
everywhere  a  maximum  has  been  reduced  to  a  quadrature.  The  assumptions  for  which 
this  solution  has  been  obtained  are  that  the  following  parameters  of  the  problem 
are  constants:  ratio  of  specific  heats  and  stagnation  conditions  of  the  gas; 
mass,  size  and  drag  coefficient  of  the  particles.  The  quadratures  (in  non- 
dimensional  form)  have  been  evaluated  for  three  values  of  the  ratio  of  specific 
heats,  namely,  5/4,  7/5,  and  5/3. 

The  results  of  this  investigation  indicate  that  a  plasma  arc  with  an  appropriately 
optimized  nozzle  appears  to  offer  the  possibility  of  significantly  Improved 
simulation  of  the  high  heating  rates  simultaneously  with  high  particle  velocities 
and  longer  testing  times  required  for  advanced  hypersonic  vehicles.  Furthermore, 
the  particle  velocities  for  a  fixed  length  of  such  an  optimized  nozzle  can  be 
Increased  by  (a)  Increasing  the  gas  stagnation  enthalpy,  the  gas  stagnation 
density,  and  the  particle  drag  coefficient;  and  (b)  decreasing  the  particle 
size  and  the  particle  density. 

Comparative  calculations,  assuming  a  particle  drag  coefficient  Cg  *  1.0,  for  17.25 
foot  long  nozzles  operating  at  2000  Btu/lbm  stagnation  enthalpy  and  75  atmospheres 
stagnation  pressure  show  that  whereas  a  conical  nozzle  may  achieve  particle 
velocities  of  6,000  and  2,400  feet  per  second  for  particle  diameters  of  50  and 
1,500  microns,  respectively,  an  optimized  nozzle  may  achieve  6,800  and  4,040 
feet  per  second  for  the  same  particle  sizes  of  50  and  1,500  microns  diameter, 
respectively.  Increasing  the  stagnation  enthalpy  to  6,000  Btu/lbm  and  the 
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stagnation  pressure  to  190  atmospheres,  the  optimized  nozzle  may  achieve  particle 
velocities  of  11,800  and  7,000  feet  per  second  for  particle  diameters  of  50 
and  1,500  microns,  respectively. 


Recommendations  for  future  work  Include  first,  the  extension  of  the  present 
analysis  to  Include  arbitrary  values  of  the  ratio  of  specific  heats  by  utilizing 
numerical  techniques  for  the  evaluation  of  the  quadrature  relating  stream 
velocity  to  the  distance  downstream  from  the  particle  Injection  point. 

Parametric  calculations  could  then  be  effeclently  performanced  to  determine 
quantitatively  the  effects  of  particle  size,  stagnation  conditions,  length  of 
device,  working  fluid  (gamma  and  molecular  weight)  on  achievable  particle 
velocities.  Next,  Investigate  numerical  techniques  to  solve  the  nozzle 
optimization  problem  with  fewer  assumptions.  In  particular,  the  assumption 
of  constant  drag  coefficient  would  be  removed  and  a  variation  of  drag  coefficient 
with,  at  least,  Reynolds  number  and  (If  sufficient  data  Is  available)  with 
Mach  number  would  be  used  In  the  determination  of  particle  forces.  Another 
approximation  that  can  be  removed  is  that  the  gas  is  calorlfically  perfect  with 
constant  y  and  constant  composition  to  Include  the  effects  of  chemical  recom¬ 
bination  and  the  variation  of  specific  heats  with  temperature.  Finally,  the 
potential  capabilities  of  existing  plasma  arc  facilities  as  particle  acceler¬ 
ators  should  be  compared.  The  characteristics  of  existing  plasma  arcs  should 
be  gathered  and  the  results  of  the  present  optimization  study,  or  the  recommended 
extensions,  should  be  applied  to  determine  the  ability  of  these  arcs  to  accel¬ 
erate  particles  when  modified  by  the  addition  of  optimum  nozzles.  The  merit  of 
this  recommended  work  would  be  either  to  make  It  possible  to  utilize  existing 
facilities  to  achieve  velocities  required  by  advanced  vehicles  by  modification 
or  extension  of  the  hardware,  or  to  lay  the  foundation  for  the  design  of  new 
hardware,  such  as  heaters  with  higher  enthalpy-pressure  capability. 
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Reduced  stream  velocity,  v„ 


Figure  1.  RELATION  BETWEEN  REDUCED  VELOCITIES  FOR  MAXIMUM  PARTICLE 
ACCELERATION 
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gure  2. 


VARIATION  OF  REDUCED  MAXIMUM  PARTICLE  ACCELERATION  WITH  REDUCED 
STREAM  VELOCITY 
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Reduced  particle  velocity. 
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Figure  10.  VARIATION  OF  PARTICLE  VELOCITY  AT  NOZZLE  EXIT  WITH  PARTICLE 
SIZE  FOR  A  17.25  FOOT  LONG  NOZZLE;  Cp  -  1.0 


